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In the framework of schematic hard spheres lattice models we discuss Edwards ' Statis- 
tical Mechanics approach to granular media. As this approach appears to hold here to 
a very good approximation, by analytical calculations of Edwards ' partition function at a 
mean field level we derive the system phase diagram and show that "jamming" corresponds 
to a phase transition from a "fluid" to a "glassy" phase, observed when crystallization is 
avoided. The nature of such a "glassy" phase turns out to be the same found in mean field 
models for glass formers. In the same context, we also briefly discuss mixing/segregation 
phenomena of binary mixtures: the presence of fluid- crystal phase transitions drives seg- 
regation as a form of phase separation and, within a given phase, gravity can also induce 
a kind of "vertical" segregation, usually not associated to phase transitions. 

1. Introduction 

We review some recent results in the study of glasses and dense granular materials 
obtained with the aid of schematic lattice models, simple enough to allow a full numerical 
and analytical investigation. Glasses and granular media exhibit deep similarities in 
their "jamming" behaviors, but although the idea of a unified description is emerging (see 
Silbert, O'Hern, Liu and Nagel contribution to this volume), the precise nature of jamming 
in non-thermal systems, such as granular media, and the origin of its close connections to 
glassy phenomena in thermal ones are still open and very important issues. 

An important conceptual open problem concerning granular media, is the absence of 
an established theoretical framework where they might be described. Edwards [ Q (see 
Edwards, Brujic, and Makse contribution to this volume) proposed a solution to such 
a problem by introducing the hypothesis that time averages of a system, exploring its 
mechanically stable states subject to some external drive (e.g., "tapping"), coincide with 
suitable ensemble averages over its "jammed states". We discuss here some recent results 
validating and generalizing Edwards' proposal (see also Dean and Lefevre, and Tarjus 
and Viot contributions). In particular, within this approach we show that mean field 
models for granular media undergo a phase transition from a (supercooled) "fluid" phase 
to a "glassy" phase, when their crystallization transition is avoided. The nature of such 
a "glassy" phase results to be the same found in mean field models for glass formers: a 
discontinuous one step Replica Symmetry Breaking phase preceded by a dynamical freez- 
ing point. These results are supported by Monte Carlo (MC) "tap dynamics" simulations 
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which, in the region of low MC shaking amphtudes, show a pronounced jamming similar 
to the one found in experiments on granular media (see Bideau, Philippe, Ribiere and 
Richard, and Caballero, Lindner, Ovarlez, Reydellet, Lanuza and Clement, and D'Anna 
and Mayor contributions to the volume). As an application of Edwards' approach to pow- 
ders we also briefly discuss segregation/mixing phenomena in these systems and discuss 
close correspondences with experiments (see Reis, Mullin and Ehrhardt contribution). 

Sec|2] gives an short introduction to the phenomenology of granular materials relevant 
to the present purpose. SecO] briefly reviews the essential lines in Edwards' Statistical 
Mechanics of powders and discuss some correspondence with other systems "frozen" at 
T = 0, such as glasses in their inherent structures. In Sections 0] to IHl we describe our 
results obtained in the study of hard spheres lattice models of granular materials, in order 
to check and apply Edwards' Statistical Mechanics to jamming and segregation. These 
are obtained by computer simulations and by analytic calculations with the generalized 
Bethe-Peierls method, or "cavity method", recently developed by Mezard and Parisi [|2]. 
Finally, we draw some conclusions. 

2. Granular Media 

Granular media are large conglomerations of discrete macroscopic particles. Ubiquitous 
in the world around us, in many respect they behave differently from any of the other 
familiar forms of matter - solids, liquids, or gases (for a review see [E]). The static 
and dynamic properties of granular media, which are strongly dissipative systems, are 
not affected by temperature, because thermal fluctuations are usually negligible. In sand 
the potential energy of a grain raised a distance equal to its own diameter can be 10^^ 
times ksT at room temperature; therefore the temperature of the external bath can be 
considered zero and these media called non-thermal. As the system cannot explore its 
phase space, unless perturbed by external forces (such as shaking or tapping), it is frozen, 
at rest, in its mechanically stable microstates. 

Even though granular media may form crystalline packings, in most cases they are found 
at rest in disordered configurations, characterized by "fluid" like distribution functions. 
Actually, gently shaken granular media exhibit a strong form of "jamming" [1110111, i-e., 
an exceedingly slow dynamics, which shows deep connections to "freezing" phenomena 
observed in many thermal systems such as glass formers [EUHl. It was discovered by the 
Chicago group [ i^, that under tapping the density of granular media tends to increase 
very slowly as a function of the number of shakes, in a stretched exponential [(HI or even 
logarithmic [ |31 way. These systems have typical "aging phenomena" along with huge 
relaxation times diverging a la Arrhenius or Vogel and Fulcher [iHiniini, similarly to glass 
formers in the freezing region. 

Here, we discuss the nature of jamming in non-thermal systems [jSUHl; and the origin 
of its close connections to glassy phenomena in thermal ones, in the framework of the 
Statistical Mechanics of powders introduced by Edwards [HHTUl: we derive the phase dia- 
gram of granular media and explain in a quantitative way their similarities and differences 
with glass formers. The introduction of a Statistical Mechanics for powders is grounded 
on observations from experiments [ (H (HI and simulations [ El that, as much as in 
thermal systems, their macroscopic properties at stationarity are characterized by a few 
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control parameters and their macrostates correspond to a huge number of microstates. 

As an apphcation of the Statistical Mechanics of powders, we also consider the intriguing 
phenomenon of segregation: in presence of shaking a granular system is not randomized, 
but its components tend to separate [ I5j. An example is the so called "Brazil nut" effect 
(BNE) where, under shaking, large particles rise to the top and small particles move 
to the bottom of the container. Interestingly, by changing grains sizes or mass ratio or 
shaking amplitudes a crossover towards a "reverse Brazil nut" effect (RBNE) was more 
recently discovered [ 120] where small particles segregates to the top and large particles 
to the bottom. Several mechanisms have been proposed to explain these phenomena 
which, although of deep practical and conceptual relevance, are still largely unknown [IT3]. 
Geometric effects, such as "percolation" [[T^ or "reorganization" [HZHIH], are known to be 
at work since, in a nutshell, small grains appear to filter beneath large ones. "Dynamical" 
effects, such as convection [ [10] or inertia [ , were shown to play a role as well. Recent 
simulations and experiments have, however, outlined that segregation phenomena can 
involve "global" mechanisms, such as "condensation" [ 120] or, more generally, "phase 
separation" [122]- We focus on these properties here. 



3. Approaches a la Edwards to Statistical Mechanics of granular media 

In the Statistical Mechanics of powders introduced by Edwards [ Hj it is postulated 
that the system at rest (i.e., not in the "fluidized" regime) can be described by suitable 
ensemble averages over its "mechanically stable" states (called here "inherent states"). 
Edwards proposed a method to individuate the probability, Pr, to find the system in its 
inherent state r, under the assumption that these mechanically stable states have the 
same a priori probability to occur. The knowledge of Pr has the conceptual advantage to 
substitute time with ensemble averages, allowing the description of the system properties 
by use of few basic theoretical concepts, as in thermodynamics. A possible approach to 
find Pr is as follows [Q] QO]- Pr is obtained as the maximum of the entropy, 

S = - J^PrlnPr (1) 

r 

with the macroscopic constraint, in the case of the canonical ensemble, that the system 
energy, E = J^r^rEr, is given. This assumption leads to the Gibbs result: 

Pr oc e-^-"-^-^'- (2) 

where Pconf is a Lagrange multiplier, called inverse configurational temperature, enforcing 
the above constraint on the energy: 

Pconf = Sconf = ^^^isi^) (3) 

Here, flis{E) is the number of inherent states with energy E. Thus, summarizing, the 
system at rest has Tbath = and Tco„/ = (3~l,f 7^ 0. 

These basic considerations, to be validated by experiments or simulations, settle a 
theoretical Statistical Mechanics framework to describe granular media. Consider, for 
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definiteness, a system of monodisperse hard spheres of mass m. In the system whole 
configuration space ^Tot, "we can write Edwards' generalized partition function as: 



where Ti,HC is the hard core interaction between grains, mgH is the gravity contribution 
to the energy {H is particles height), and the factor 11,. is a projector on the space of 
"mechanically stable" states (i.e., inherent states space) flis- if r G flis then 11^ = 1 else 
Ur = 0. Usual Statistical Mechanics, where is identified with Tbath, is recovered in 
the case where 11^ = 1 Vr. 

In order to test Edwards proposal one introduces a dynamics (such as "trapping" ) that 
allows the system to explore the inherent states space. Using this djTiamics one has first 
to check that at stationarity the system properties do not depend on the details of the 
dynamical history, i.e., a "thermodynamic" description is indeed possible. Then one must 
check that time averages obtained using such a dynamics compare well with ensemble 
averages over the distribution Eq.Q. 

3.1. Configurational entropy of granular media and glasses frozen at T = 

Edwards approach could also be relevant for a supercooled liquid quenched at very 
low temperature (about zero). In this case the system remains blocked in its "inherent 
structures", i.e., the local minima of the potential energy in the particle configuration 
space [ 1221 1211 UH 1121 (by analogy with the glass terminology we call "inherent states" 
the mechanically stable states in granular materials). 

In glasses, in the inherent structure approach [ (23 1211 ^1 the configurational 
entropy (and consequently the configurational temperature) associated to the number of 
inherent structures corresponding to a given energy, E, can be defined too. In this context 
two different ways were essentially used to allow the system explore its inherent structure 
space. One way is by quenching the system over and over from an equilibrium tempera- 
ture, T, to zero temperature [I231IIIIII21- Using this procedure, Sciortino et al. [[TT] found 
that the configurational temperature numerically coincides with the equilibrium tempera- 
ture T, provided that T is low enough. Another way to visit the inherent structures is by 
letting the system aging in contact with an almost zero bath temperature, Ti,ath- During 
the aging process an effective temperature, T^yn, can be defined via the off-equilibrium 
extension of the Fluctuation-Dissipation ratio [I2ni- It happens that in mean field models 
[125, J6J this effective temperature coincides with the configurational temperature, Tconf- 
This was also shown in Ref. [l2Hj5 in the limit Ti,ath — > 0, for a class of finite dimensional 
systems. The connection between Edwards approach for granular media and the results 
in glass theory has been pointed out in [|271I2H1 Cll 1221 HHj- 

A second procedure, called "tap dynamics", similar to that used in the compaction of 
real granular materials was instead used in [ 123 HDl 1201 1211 1221 1221 H] • In particular in 
Ref. [ 123 ^1 each tap consists in raising the bath temperature to a value Tr and, after 
a lapse of time tq, quenching it back to zero. By cyclically repeating the process the 
system explores the space of the inherent states. Once the stationary state is reached a 
temperature, T/^, is defined via the equilibrium Fluctuation-Dissipation relation, and if 
Edwards' assumption applies, Tfd coincides with the configurational temperature. This 
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Figure 1. Our lattice models for granular media are subject to a Monte Carlo dynamics made of "taps" 
sequences. A "tap" is a period of time, of length tq (the tap duration), during which the system evolves 
at a finite bath temperature Tr (the tap amplitude); after each "tap" the system evolves at Tr = and 
reaches an inherent state. 



was verified for different finite dimensional models [123 CD]- Many other studies discussing 
Edwards' approach in this perspective were also presented [ EOl EH E2] • In the following 
sections we give a picture of results [1211110] about such a statistical mechanics approach 
in schematic models of granular media and glasses. 

4. Lattice models for granular media 

4.1. A monodisperse hard-sphere system under gravity 

The simplest model for granular media we considered [ 10 is a system of hard-spheres 
of equal diameter = 1, subjected to gravity. We have studied this model on a lattice, 
constraining the centers of mass of the spheres on the sites of a cubic lattice (see inset in 
Fig. E)). The Hamiltonian of the system is: 

H = nHci{n'i}) + gm'^niZi, (5) 

i 

where the height of site i is Zi, g = 1 is gravity acceleration, m = 1 the grains mass, 
nj = 0, 1 the usual occupancy variable (i.e., nj = or 1 if site i is empty of filled by 
a grain) and ?i//c({^i}) an hard-core interaction term that prevents the overlapping of 
nearest neighbor grains (this term can be written as Ti-Hciini}) = JJ2{ij) ^i^jy where the 
limit J ^ oo is taken). 

The grains are subject to a dynamics made of a sequence of Monte Carlo "taps" (see Fig. 
P): a single "tap" is a period of time, of length tq (the tap duration), where particles can 
diffuse laterally, upwards with probability pup G [0, 1/2], and downwards with probability 
1 — Pup- When the "tap" is off grains can only move downwards (i.e., pup = 0) and the 
system evolves with pup = until it reaches a blocked configuration (i.e., an "inherent 
state") where no grain can move downwards without violating the hard core repulsion. 
The parameter pup has an effect equivalent to keep the system in contact (for a time tq) 
with a bath temperature Tr = mgao/ \n[{l — Pup)/Pup] (called the "tap amplitude"). The 
properties of the system are measured when this is in a blocked state. Time averages. 
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Figure 2. The time average of the energy, E, 
and (inset) its fluctuations, /S.E , recorded at 
stationarity during a tap dynamics, as a func- 
tion of the tap amplitude, Tr, in the 3D lat- 
tice monodisperse hard sphere model. Different 
curves correspond to sequences of tap with dif- 
ferent values of the duration of each single tap. 
To- 
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Figure 3. Time averages of energy fluctua- 
tions Ai?2 plotted as function of the time aver- 
age of energy E. k and ■ are time averages 
obtained with different tap dynamics. ^'^^ 
independently calculated ensemble averages ac- 
cording to Eq.(|2Jl. The collapse of the data ob- 
tained with different dynamics shows that the 
system stationary states are characterized by a 
single thermodynamic parameter. The agree- 
ment with the ensemble averages show the suc- 
cess of Edwards' approach to describe the sys- 
tem macroscopic properties. 



therefore, are averages over the blocked configurations reached with this dynamics. Time 
t is measured as the number of taps apphed to the system. 

Under such a tap dynamics the systems reaches a stationary state where the Statistical 
Mechanics approach to granular media can be tested, and particularly Edwards hypothesis 
can be verified by comparing time averages to ensemble averages of Eq.®. 

4.2. The stationary states of the tap dynamics 

During the tap dynamics, in the stationary state, the time average of the energy, E, and 
its fluctuations, AE"^, are calculated. Figure El shows E (main frame) and /S.E (inset) 
as function the tap amplitude, Tr, (for several values of the tap duration, tq). Since 

2 

sequences of taps, with same Tr and different tq, give different values of E and Ai? , it 
is apparent that Tp is not the right thermodynamic parameter. On the other hand, if 
the stationary states are indeed characterized by a single thermodynamic parameter the 
curves corresponding to different tap sequences (i.e. different Tp and tq) should collapse 

2 

onto a single master function, when AE is parametrically plotted as function of E. This 
is the case in the present model, where the data collapse is in fact found and shown in 
Fig. El This is a prediction that could be easily checked in real granular materials. 
A technique to derive from raw data the thermodynamic parameter (3fd conjugated to 
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Figure 4. Main frame The time average E 
and the ensemble average over the distribution 
Eq.jSJ {E), plotted respectively as a function 
of Tfd and Tconf (in units mgao), in the 3D 
monodisperse hard-sphere system under grav- 
ity described in the text. Symbols are as in Fig. 
121 Time averages over the tap dynamics and 
Edwards' ensemble averages coincide. Lower 
Inset The temperature Tfd = Pjd defined by 
Eq.JHJ as function of Tr (in units mgao) for 
To = 500,10,5 MCS (from top to bottom). 
The straight line is the function Tfd = Tr. 



Figure 5. The density self-overlap function, 
Q, and (upper inset) the system density on 
the bottom layer, pb, plotted as function of 
Tfd (in units mgao), compared with the en- 
semble averages over the distribution Eq.© 
(the black empty circles), plotted as function 
of Tconf (in units mgao), in the 3D monodis- 
perse hard-sphere system under gravity. Also 
for this system there is a very good agreement 
between the independently calculated time av- 
erages over the tap dynamics and the statistical 
mechanics ensemble averages a la Edwards. 



E (apart from an integration constant, (3q), is through the usual equilibrium Fluctuation- 
Dissipation relation: 

dE 



AE^. (6) 



By integrating Eq.®, Pfd — Po can be expressed as function of E or (for a fixed value of 
To) as function of /5r = l/^r: Pfd — Po = diPr) (the constant /3o can be determined as 
explained in [CD])- 

4.3. Edwards' averages 

In Sect. 14.21 we have found that the fluctuations of the energy in the stationary state 
depend only on the energy, E, and not on the past history. More generally, we found 
[ Uni that all the macroscopic quantities we observed depend only on the energy, E, 
or on its conjugate thermodynamic parameter, Pf^, thus the stationary state can be 
genuinely considered a "thermodynamic state" . In this case one can attempt to construct 
an equilibrium statistical mechanics, as originally suggested by Edwards [Hj. 

To test Edward hypothesis one has to compare the time average of the energy, E{i3fd), 
recorded during the taps sequences, with the ensemble average, {E){Pconf), over the dis- 
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tribution Eq.©. Therefore we have independently calculated the ensemble average {E), 
as function of Pconf- Fig. HI (see also Fig. |2)) shows a very good agreement between 
{E){Pconf) and E{l3fd) (notice that there are no adjustable parameters). In order to check 
the scenario in further details, we have also calculated the system density on the bottom 
layer, p^, and the density self-overlap function Q: Q = ^/Nj^iLi ^f^^f \ where the sum 
is over all the lattice sites and n^^-* and n^"^^ refer to two different copies of the systems, 
which are at stationarity under the same tap dynamics. We have verified that, when 
plotted as function of Tj^, both pi, and Q scale on a single master function (see Fig. E}. 
In Fig. m (inset) we also show the dependence of the configurational temperature Tconf 
on the parameters of the tap dynamics Tr and tq. Finally, we mention that we have 
also successfully tested Edwards scenario in an other model, the "frustrated lattice gas" 
[ [TU] , a system in the category of spin glasses. 



4.4. Relaxation during MC "tap" dynamics 

The MC tap dynamics, in both the "frustrated lattice gas" and hard sphere models, 
exhibits a rich structure in agreement with experimental findings [I^IHl- In the region of 
small tap amplitudes, the system gets "jammed" and "memory" phenomena are observed, 
along with "irreversibility" effects [ El El- In particular, it is interesting to consider 
correlation functions such as C(t,tw) = B{t,ty^) / Bitw^t^), where B(t,tw) = + 
tw)ni{tw)) — {riiit + tw)) {riiitw))]. In the high Tr region, C{t,tyj) has a time translation 
invariant (TTI) behavior, i.e., C{t,tw) = C{t). Asymptotically C{t) can be well fitted by 
stretched exponentials: C{t) = Coexp[— (t/r)^]. The exponent /3 becomes significantly 
lower than 1 at low amplitudes. The above fit defines the relaxation time r(Tr):the growth 
of r by decreasing Tr is well approximated by an Arrhenius or Vogel-Tamman-Fulcher 
law (as early found in [EHIS]), resembling the slowing down of glass formers close to the 
glass transition, a result also recently experimentally reported in granular media [Elini- 
r ~ To exp [Eo/{Tr — T^)] ■ The divergence point, (which in simulations is difficult to 
precisely locate; here it is between and 1), of r is interpreted as the numerical location of 
the point of dynamical arrest of the system, where an "ideal" transition to a glassy phase 
occurs. By quenching the system at low values of Tp, the TTI character of relaxation is 
lost and logarithmic aging behaviors are found. For slow quenches the hard spheres model 
is able, anyway, to attain its crystal phase. The precise nature of the "glassy" region is, 
however, very difficult to be numerically determined, so in the next section an analytic 
approach is presented. 

An other dynamical quantity of interest is grains mean square displacement, R'^{t), 
defining for t —>■ oo the diffusion coefficient D{p) = R^{t)/t. As much as in glass formers, D 
goes to zero around a high density pc, signaling a localization transition in which particles 
are confined in local cages and the macroscopic diffusion-like processes are suppressed [ 
|7]. This phenomenon may also be described in a different way: pc is the density above 
which it becomes impossible to obtain a macroscopic rearrangement of particles positions 
without increasing the system volume, i.e., the density at which macroscopic shear in the 
system is impossible without dilatancy. This remark outlines the close correspondence 
with the phenomenon of the Reynolds transition in granular media [ [7] . 
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5. A mean field study of hard spheres under gravity on a random graph 

We have discussed up to here schematic models for granular media shown to be well 
described by Edwards' assumption. Since the exact calculation of Edwards partition 
function, Z, for the above model of monodisperse hard spheres on a cubic lattice is hardly 
feasible, we now discuss a mean field theory [EH]- In an approximation a la Bethe- 
Peierls, we consider a random graph version of such a lattice, sketched in Fig. IHl more 
specifically we consider a 3D lattice box with H horizontal layers (i.e., z G 
occupied by hard spheres; each layer is a random graph of given connectivity, k — 1 (we 
take /c = 4); each site in layer z is also connected to its homologous site in z — 1 and 
z + 1 (the total connectivity is thus + 1). The Hamiltonian is the one of Eq.© plus 
a chemical potential term to control the overall density. Hard Core repulsion prevents 
two connected sites to be occupied at the same time. In the present lattice model we 
adopt a simple definition of "mechanical stability": a grain is "stable" if it has a grain 
underneath. The operator 11^ has thus a simple expression: 11^ = lim^^oo exp {—KTiEdw} 
where HEdw = Y.i ^n,{z),iSn,{z-i),oSn,{z-2)fi (for clarity, we have shown the z dependence 
in ni{z)). 




z+l 



z-1 



Figure 6. In our mean field approximation, the hard spheres describing the system grains are located on 
a Bethe lattice, sketched in the figure, where each horizontal layer is a random graph of given connectivity. 
Homologous sites on neighboring layers are also linked and the overall connectivity, c, of the vertices is 
c = fc + 1 = 5. 

By using the "cavity method" [i2j, the phase diagram is found [ESJ- At low Ns or 
high Tconf a fiuid-like phase is found, characterized by a homogeneous Replica Symmetric 
(RS) solution, in which only one pure state exists and the local fields are the same for 
all the sites of the lattice (translational invariance). For a given Ng, by lowering T^onf 
(see Figs. [7|and|H)), a phase transition to a crystal phase (an RS solution with no space 
translation invariance) is found at Tm- Notice that the fiuid phase still exists below Tm 
as a metastable phase corresponding to a supercooled fiuid found when crystallization 
is avoided. Within the one-step replica symmetry breaking (IRSB) ansatz of the cavity 
method [ Ej, a non trivial solution appears for the first time at a given temperature 
Td{Ns), signaling the existence of an exponentially high number of pure states. In mean 
field theory To is interpreted as the location of a purely dynamical transition as in mode- 
coupling theory, but in real systems it might correspond just to a crossover in the dynamics 
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Figure 7. The system mean field phase di- 
agram is plotted in the plane of its two con- 
trol parameters {Tconf, Ng): T^onf is Edwards' 
"configurational temperature" and Ns the av- 
erage number of grains per unit surface in the 
box. At low Ng or high Tconf, the system 
is found in a fluid phase. The fluid forms a 
crystal below a melting transition line T^iNs). 
When crystallization is avoided, the "super- 
cooled" (i.e., metastable) fluid has a thermo- 
dynamic phase transition, at a point Tk{Ns), 
to a Replica Symmetry Breaking "glassy" phase 
with the same structure found in mean field the- 
ory of glass formers. In between Tm(iVs) and 
TKiNs) a dynamical freezing point, ToiNg), is 
located, where the system characteristic time 
scales diverge. 
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Figure 8. For a system with a given number 
of grains (i.e., a given Ng), the overall num- 
ber density, <& = Ng/2{z) {{z) is the average 
height), calculated in mean field approximation 
is plotted as a function of Tconf, ^(Tconf) has 
a shape very similar to the one observed in the 
"reversible regime" of tap experiments and MC 
simulations of the cubic lattice model for <i>(rr). 
The location of the glass transition, Tk (filled 
circle), corresponds to a cusp in the function 
^{Tconf)- The passage from the fluid to su- 
percooled fluid is Tm (filled square). The dy- 
namical crossover point Td is found around the 
flex of ^ {Tconf) and well corresponds to the po- 
sition of a characteristic shaking amplitude F* 
found in experiments and simulations where the 
"irreversible" and "reversible" regimes approx- 
imately meet. 



(see [ EEl Cni inni and Ref.s therein). The IRSB solution becomes stable at a lower point 
Tk, where a thermodynamic transition from the supercooled fluid to a IRSB glassy phase 
takes place (see Fig. [7j) in a scenario a la Kauzmann with a vanishing complexity of pure 
states (which stays flnite for Tk <T < Td). 

The results of these calculations, summarized in the phase diagram of Fig. [7[ are further 
illustrated in Fig. |H1 in a system with a given number of grains (i.e., a given A^^), the 
overall number density, $, is plotted as a function of Tconf (here by deflnition $ = Ns/2{z), 
where (z) is the average height). The shown curve, ^{Tconf), is the equilibrium function 
here calculated. It has a shape very similar to the one observed in tap experiments [ 
miE], or in MC simulations on the cubic lattice (see also [Ej), where the density is plotted 
as a function of the shaking amplitude F (along the so called "reversible branch"). In 
particular, a comparison of our mean fleld results with simulations of the 3D model of 
Hard Spheres under the tap dynamics shows a very good agreement. 
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Summarizing, in the present mean field scenario of a granular medium with Ns particles 
per surface, in general, at high Tconf (i-e. high shaking amplitudes) a fluid phase is 
located (see Fig. [7j). By lowering Tconf, a phase transition to a crystal phase is found 
at Tm- However, when crystallization is avoided, the fluid phase still exists below 
as a metastable phase corresponding to a supercooled fluid. At a lower point, Tq, an 
exponentially high number of new metastable states appears, interpreted, at a mean field 
level, as the location of a purely dynamical transition, which in real system is thought 
to correspond just to a dynamical crossover. Finally, at a even lower point, T^, the 
supercooled fluid has a genuinely thermodynamics discontinuous phase transition to glassy 
state. The structure of the glass transition of the present model for granular media, 
obtained in the framework of Edwards' theory, is the same found in the glass transition 
of the p-spin glass and in other mean field models for glass formers [ 1^ I13j. 

6. A hard-sphere binary mixture under gravity 

In order to study segregation and to test Edwards proposal to a more complicate system, 
here we consider a hard-sphere binary system made of two species 1 (small) and 2 (large) 
with grain diameters oq and \/2ao, under gravity on a cubic lattice of spacing oq = 1. We 
set the units such that the two kinds of grain have masses mi = 1 and m2 = 2mi, and 
gravity acceleration is g = 1. The hard core potential Hhc is such that two large nearest 
neighbor particles cannot overlap. This implies that only couples of small particles can 
be nearest neighbors on the lattice. The system overall Hamiltonian is: 

^ = 'Hhc + migHi + m2gH2, (7) 

where Hi = Zi and H2 = Ylf'^ height of site i is Zi and the two sums are over 

all particles of species 1 and 2, respectively. In the above units, the gravitational energies 
in a given configuration are thus Ei = Hi and E2 = 2H2. 

As before, grains are confined in a box of linear size L between with periodic boundary 
conditions in the horizontal directions and initially prepared in a random loose stable 
pack. Under the tap dynamics the system approaches a stationary state for each value of 
the tap parameters Tr and tq used. In Fig. "we plot as function of Tr (for several values 
of Tq) the asymptotic value of the fertica/ segregation parameter, i.e., the difference of the 
average heights of the small and large grains at stationarity, A/;,(Tr,ro) = hi — h2- Here 
hi and h2 are the average of Hi/Ni and H2/N2 over the tap dynamics at stationarity. 

The results given in the main panel of Fig. IHl apparently show that Tr is not a right 
thermodynamic parameter, since sequences of taps with different tq give different values 
for the system observables. However, if the stationary states corresponding to different tap 
dynamics (i.e., different Tp and Tq) are indeed characterized by a single thermodynamic 
parameter, as in the monodisperse case above, the curves of Fig. El should collapse onto 
a universal master function when A/i(Tr,ro) is parametrically plotted as function of an 
other macroscopic observable such as the average energy, e(Tr,ro) = {Ei + E2)/N {N is 
the total number of particles). This collapse of data is not observed here, as it is apparent 
in the inset of Fig. El We found, instead, [ TI? that two macroscopic quantities can be 
sufficient to characterize uniquely the stationary state of the system. These two quantities 
are, for instance, the energy e and the height difference Ah. Of course since e = ahi + 2bh2 
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Figure 9. Main frame The difference of 
the average heights of small and large grains, 
Ah = hi — h2, measured at stationarity in the 
binary hard spheres mixture under gravity, is 
plotted as a function of tap amplitude, Tr (in 
units mgao). The three sets of points corre- 
spond to the shown tap durations, tq. At high 
Tr larger grains are found above the smaller, 
i.e. Ah < 0, as in the Brazil nut effect (BNE). 
Below a T-^{to) the opposite is found (Reverse 
Brazil nut effect, RBNE). Inset The Ah data 
of the main frame are plotted as a function of 
the corresponding average energy, e. The three 
sets of data do not collapse, as before, onto a 
single master function. 



Figure 10. Main frame The average den- 
sity of large grains on the box bottom layer, 
measured at stationarity for different Tr and tq, 
scale almost on a single master function when 
plotted as a function of the large grains height, 
/i2- Upper inset The average number of con- 
tacts between large grains per particle, Nc, ob- 
tained for different Tr and tq, scale on a single 
master function when plotted as a function of 
the system energy, e. 



(where a = Ni/N and b = N2/N) and Ah = /ii — /12, we can also choose hi and /i2 to 
characterize the stationary state. Namely, we found that a generic macroscopic quantity 
A, averaged over the tap dynamics in the stationary state, is only dependent on hi and 
/12, i-e., A = A{hi,h2)- We have checked that this is the case for several independent 
observables, such as the number of contacts between large particles, Nc, the density of 
small and large particles on the bottom layer, p\ and and others. In particular, as 
shown in Fig. we find with good approximation that: Nc — Nc{e) = Nc{ahi + bh2), 
p\ ^ P2(^2)) Pi — Pi(^i)- Therefore we need both hi and h2 to characterize unambiguously 
the state of the system; namely all the observables assume the same values in a stationary 
state characterized by the same values of hi and /12, independently on the previous history 
(i.e., in our case independently on the particular tapping parameters Tr and tq). 

These findings imply that an extension of Edwards' approach is required, where at least 
tvjo thermodynamic parameters have to be included [Cni- As before, this can be obtained 
by assuming that the microscopic distribution is given by the principle of maximum 
entropy with the constraint that the average gravitational energies of the two species 
El = J2r PrEir and E2 = PrE2r are independently fixed. This gives two Lagrange 
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multipliers: 

^ _ dlnQis{Ei,E2) . _ d\nQjs{E^,E2) 

— dE[ — m 

where Qjs{Ei, E2) is the number of inherent states with Ei,E2. 

The hypothesis that the ensemble distribution at stationarity is the above can be tested 
as we have already previously shown. We have to check that the time average of any quan- 
tity, A{hi, h2), as recorded during the taps sequences in a stationary state characterized 
by given values hi and /12, coincides with the ensemble average, {A){hi, h2), over the 
generalized version of distribution Eq.Q. To this aim, we have calculated the ensemble 
averages (N^), {p\), {p\) for different values of /?i and P2', we have expressed parametrically 
(iVc), (^2); (Pi); as function of the average of hi and h2, and compared them with the 
corresponding quantities, Nc, p\ and averaged over the tap dynamics. The two sets 
of data are plotted in Fig. [TUl showing a good agreement (notice, there are no adjustable 
parameters) . 
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Figure 11. Mean field phase diagram of the binary mixture in the plane (iVi^, N2s) (particles densities 
per unit surface of species 1 and 2) for mi(3i = 0.8 and m2/32 — 1-25 (main frame). A fluid and a 
crystal phase are present (glassy and metastable phases are not shown here). Size segregation is found 
as a form of phase separation when the system is prepared out of these phases. Within a given phase, 
gravity drives vertical mixing/segregation phenomena (the BNE, RBNE and MIXING regions) separated 
by smooth crossovers (shaded areas). Inset Phase diagram for miPi — 1.25 and TO2/32 — 0.8. Here also a 
new phase, with strong BNE, appears where a fluid, rich in small grains, is found beneath a crystal rich 
in large ones ("C+F" region). For mi/3i << TO2/32, strong RBNE can be found, with the fluid above the 
crystal. 

In [I2Z] the mean field approximation of Sec. El is applied to the present binary mixture, 
giving a precise interpretation of segregation phenomena observed in the model of Eq.((7|). 
as shown in FigHH In a nutshell, by changing A^'i^, N2S (i.e., densities per unit surface 
of species 1 and 2), j3i and (32 (i.e, their configurational temperatures), or masses or 
sizes ratios, the system exhibits true phase transitions from fluid to crystal phases. As 
much as in thermal media, this induces segregation effects associated to phase separation 
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phenomena, with the formation of coexisting phases rich in small or large grains, as 
experimentally observed in [22J. Within a given phase, the presence of gravity can also 
drive a form of "vertical" segregation, where large grains are found on average above (the 
well known BNE) or below small grains (RBNE [,20^), but this is usually not directly 
associated to phase transitions. 

7. Conclusions 

Summarizing, within the schematic framework of lattice models, we have shown, by 
Monte Carlo simulations of "tap" dynamics, that Edwards' approach to granular media 
at rest appears to be well grounded. The system stationary states are indeed independent 
on the sample history as in a "thermodynamics" system, and can be described in terms of 
a distribution function characterized by a few control parameters (such as configurational 
temperatures). By use of Edwards approach, we have derived, by analytical calculations 
at a mean field level, the phase diagram of these systems. In particular, we discovered 
that "jamming" corresponds to a phase transition from a "fluid" to a "glassy" phase, 
observed when crystallization is avoided. Interestingly, the nature of such a "glassy" 
phase turns out to be the same found in mean field models for glass formers. In the same 
framework, we have also briefly discussed segregation patterns observed in a hard sphere 
binary model under gravity subject to sequences of taps. Here, the presence of fluid-crystal 
phase transitions in the system drives segregation as a form of phase separation. Within 
a given phase, gravity can also induce a kind of "vertical" segregation, not associated to 
phase transitions. 

In practice, even though the general validity of Edwards approach to "frozen" systems 
has just begun to be assessed, it turns out that a first reference framework is emerg- 
ing to understand the physics of granular media and their deep connections with other 
"jamming" systems such as glass formers. 
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